A potential WIMP signature for the caustic ring halo model 
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Weakly Interacting Massive Particle (WIMP) direct detection event rate calculations usually rely 
on fairly simple, essentially static, analytic halo models. This is largely since the resolution of 
numerical simulations is not yet large enough to allow the full numerical calculation of the WIMP 
density and velocity distribution. In this paper we study the direct detection rate, in particular its 
energy dependence and annual modulation, for the caustic ring halo model. In this model, which 
uses simple assumptions to model the infall of dark matter onto the halo, the distribution of the 
cold dark matter particles at the Earth's location has a series of peaks in velocity space. We find 
that the recoil energy spectrum contains distinctive steps and the sign of the annual modulation in 
the event rate changes as a function of recoil energy. These effects provide a potentially distinctive 
experimental signal. 
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I. INTRODUCTION 

Weakly Interacting Massive Particle (WIMP) direct 
detection experiments are just reaching the sensitivity 
required to probe the interesting range of mass-cross- 
section parameter space where relic neutralinos could 
constitute the dark matter. The DAMA collaboration, 
using a detector consisting of radiopure Nal crystal scin- 
tillators at the Gran Sasso Laboratory, have reported the 
detection of a 4<r annual modulation signal in their di- 
rect detection experiment, consistent with WIMP scat- 
tering J0,|). Whilst this result is somewhat controver- 
sial fffljfy it illustrates the potential of current and up- 
coming WIMP direct detection experiments. 

Event rate calculations and detection strategies for 
particle physics dark matter candidates are usually 
based on the assumption of a standard Maxwellian halo 
model J|-01 ■ The standard Maxwellian halo model has a 
number of deficiencies, in particular the real halo con- 
tains substructure and is not perfectly spherical and 
isotropic Since the resolution of numerical simula- 
tions is not yet large enough to allow the full numerical 
calculation of the WIMP density and velocity distribu- 
tion, analytic, or at least semi-analytic, models for the 
dark matter halo must be used. The direct detection 
rate and in particular its annual modulation, which oc- 
curs due to the Earth's motion, has been calculated for a 
range of analytic non-standard halo models [5|j9|-[l5|| . It 
has been found that the region in the mass-cross section 
plane selected by the DAMA data depends substantially 
on the halo model assumed |I(||l2],[l5| . 

Analytic halo models usually assume an essentially 
static halo, whereas in reality the halo is forming via the 
ongoing infall of surrounding dark matter |l6| ] . The caus- 
tic ring halo model, which arises from simple assumptions 
about the infall of dark matter onto the halo, provides an 
analytic model of some features of the dark matter dis- 



tribution which may result from this accretion process. 
It is therefore worthwhile to calculate the observational 
features which the caustic ring halo model produces. The 
directional WIMP direct detection rate, which would be 
probed by proposed experiments such as DRIFT |jl7) , has 
been calculated by Copi, Han and Krauss JTq] , whilst 
Vergados |l9| has calculated the total WIMP direct de- 
tection rate. In this paper we study the variation of the 
differential direct detection rate, in particular its annual 
modulation, with detector recoil energy. 



II. CAUSTIC RING HALO MODEL 

Cold dark matter (CDM) particles are collisionless and 
have low velocity dispersion (< 30kms _1 ) so that parti- 
cles falling onto an isolated galaxy are expected to oscil- 
late in and out of the galaxy a number of times before 
they are virialised by inhomogeneities (such as molecu- 
lar clouds, globular clusters and stars) |20|. These non- 
virialiscd CDM flows lead to the formation of caustic 
rings at the points where the particles with the most 
angular momentum in a given inflow reach their point of 
closest approach to the galactic centre, hence the name 
of the model. Furthermore the distribution of the parti- 
cles at any given location is expected to have a series of 
peaks in velocity space, corresponding to particles which 
are falling into the galaxy for the first time and those 
which have fallen in and out a number of times but have 
not yet been thermalized. Whilst this model is obviously 
a simplification of the hierarchical accretion process via 
which the galactic halo forms, in particular the Milky 
Way is not an isolated galaxy, the resolution of N-body 
simulations is not yet large enough to resolve these sorts 
of features. 

The velocities and densities at the Earth's location ex- 
pected due to these flows have been calculated using the 
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TABLE I. The density and velocity components, in the 
rest frame of the galaxy, of the velocity flows. 



self-similar infall model |21[ generalised to take into ac- 
count the angular momentum of the CDM particles [^2| . 
The Earth is located between the 4th and 5th caustic 
rings and the velocity flows corresponding to these two 
rings constitute roughly 30% of the local halo density. 
Analysis of 32 extended galactic rotation curves has pro- 
vided some evidence for the 1st and 2nd caustic rings [f23| , 
whilst analysis of an IRAS map of the galactic disk ap- 
parently reveals the presence of the 5th ring [^4| . 

The velocity distribution function of the velocity flows 
can be written as: 



/(v) = ZjpjSiv- Vj) 



(1) 



where pj and Vj are the density and velocity of the j-th 
flow. Table I contains the most recently calculated values 
of pj and v 3 p5| (note that there are two, inward and 
outward, flows for each velocity peak). The total density 
is po = 102gcm~ 3 = 0.57GeVcm -3 , with the velocity 
flows contributing 65% of the total. We will assume that 
the thermalized background distribution is a Maxwellian 
with velocity dispersion vq = 220 kms -1 . 



III. ANNUAL MODULATION SIGNAL 

The WIMP detection rate depends on the speed dis- 
tribution of the WIMPs in the rest frame of the detector, 
/„. This is found from the halo velocity distribution, / (v) 
by making a Galilean transformation v — * v = v — v , 
where v e is the Earth's velocity relative to the galactic 



rest frame, and then integrating over the angular distri- 
bution. In galactic co-ordinates the axis of the ecliptic 
lies very close to the (f> — z plane and is inclined at an 
angle 7 « 29.80° to the <f> — r plane. Including all com- 
ponents of the Earth's motion, not just that parallel to 
the galactic rotation jll| : 

v = v\ sin a f + 

(vq + v\ cos a sin 7) <f> — v± cos a cos 7 z , (2) 

where vo 1=3 232kms _1 is the speed of the sun with respect 
to the galactic rest frame, vi ss 30kms _1 is the orbital 
speed of the Earth around the Sun and a = 2ir(t — to)/T, 
with T = 1 year and to ~ 153 days (June 2nd), when the 
component of the Earth's velocity parallel to the Sun's 
motion is largest. 

In the range of masses and interaction cross sections 
accessible to current direct detection experiments the 
best motivated WIMP candidate is the neutralino, for 
which the event rate is dominated by the scalar contri- 
bution. The differential event rate simplifies to (see e.g. 
Refs. fUTsI for details): 



dR „ 
dE = ^ 



P0.3 (m p + m x ) 



'TTVq 



xJ A 2 T(E)F 2 (E) 



(3) 



where it is conventional to normalise the local WIMP 
density p x , to a fiducial value po.3 — 0.3GeVcm -3 , such 
that £ = p x //Oo.3: E is the energy deposited in the detec- 
tor, A is the atomic number of the detector nuclei, F(E) 
is the detector form factor (the Saxon Woods form factor 
is used for / whilst that of Na is taken to be unity, see 
e.g. Ref. @) and T(E) is defined as § 



T(E) 



/7TUo 



dv . 



(4) 



where v n 



is the minimum detectable WIMP velocity 

1/2 



E(m x + rriA)' 
2m x mA 



(5) 



m x is the WIMP mass and wla is the atomic mass of the 
target nuclei. 

In order to compare the theoretical signal with that 
observed we need to take into account the response of the 
detector. The electron equivalent energy, E cc , which is 
actually measured is a fixed fraction of the recoil energy: 
-Ecc = IkE. The quenching factors for I and Na are 
<7i = 0.09 and <7Na = 0.30 respectively p6[ . The energy 
resolution of the detector (9) is already taken into account 
in the data released by the DAMA collaboration. 

The expected experimental spectrum per energy bin 
for the DAMA collaboration set-up is then given by |T(]] 
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TABLE II. The density, (f> velocity component in the rest 
frame of the galaxy v$ , <j> velocity component in the rest frame 
of the Earth v^, and total speed in the rest frame of the 
Earth vtot, of the caustic flows in June, when a = (and in 
December when a = n). 



where rN a = 0.153 and n = 0.847 are the mass fractions 
of Na and I respectively. Since Vq 3> «i the differential 
event rate in the k-th energy bin can be expanded in a 
Taylor series in cos a B : 



Ai? 

— (E k ) ps 5 , k + ^m, k cos a . 



IV. RESULTS 



(7) 



Whilst all 3 components of the Earth's velocity need 
to be included to calculate the annual modulation sig- 
nal accurately, the signal is largely determined by the 
component in the galactic plane || : 



[1.05 + 0.06 cos a] 



(8) 



where u c ; rc = 220 kms -1 is the local circular velocity 
about the galactic centre. Before presenting the results of 
a numerical calculation, using all three components of the 
Earths' velocity, we will carry out a simple analytic cal- 
culation, using only the component in the galactic plane, 
in order to elucidate the physical origin of the variation 
in T{E). 

In June when a = 
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FIG. 1. The value of T(E) in June (solid line) and De- 
cember (dotted line) due to the velocity flows alone for four 
values of the WIMP mass m x = 30, 50, 100, 200 GeV, (top 
left, bottom left, top right and bottom right respectively). 
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FIG. 2. The value of T(E) in June (solid line) and De- 
cember (dotted line) for a halo model with caustics plus an 
isothermal background, for four values of the WIMP mass 
m x = 30, 50, 100, 200 GeV, (as before). 



whilst in December when a = it 



u c ,4> 



= 0.99 xv s = 217.8kms~ 



(10) 



w c .0 = 1.11 x v Q = 244.2kms 



(9) 



Table 2 contains the density, 0- velocity component, in 
the rest frames of the galaxy and Earth, and the total 
velocity in the rest frame of the Earth of the velocity flows 
for a — and it. In both cases the total density in flows 
with negative v$ (incident from the forward direction) 
is 17.1 xl0 -26 gcm~ 3 whilst the total density in flows 
with positive v<f> (incident from the backward direction) 
is 49.4 x 10 - 26 gem -3 i.e. there are more WIMPs incident 
from backwards than forwards as found by Copi, Han and 
Krauss jL8| . This is the opposite of the directional signal 
produced by a pure Maxwellian halo. 

In order to illustrate how the variations due to the 
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caustics are smoothed out by the isothermal background 
we plot T(E), as a function of E, for a Ge 76 detector in 
Fig. 1 for the velocity flows alone, and in Fig. 2 for the 
complete halo model described above in Sec. ||, where the 
velocity flows contribute 65% of the local density with the 
remaining 35% in an isothermal background. Values for 
other monatomic detectors can be found by rescaling the 
x-axis by toa/ (wa + m x) 2 - 

For a pure Maxwellian halo the signal is largest in 
December for small recoil energies, switching to become 
largest in June as the recoil energy is increased |ll| . The 
signal for the velocity flows alone is more complicated. 
The contribution of the j-th velocity flow to T(E) is pro- 
portional to Pj/vt t if Vtot > Vmin and is zero otherwise. 
The contribution of the high density flows to the signal 
is largest in June, since their w to t is smaller in June than 
in December. Therefore at low energies, where all the ve- 
locity flows contribute to the signal, the signal is largest 
in June. A given high density velocity flow stops con- 
tributing to the signal, vtot < v min , for smaller v m i n , or 
equivalently E, in June compared to December however. 
This means that the steplike decreases in T(E), which 
arise when a given flow stops contributing to the signal, 
occur at lower energies in June than in September (see 
Fig. 1). In other words for some range of recoil energies 
a given flow contributes to the signal in December but 
not in June. At large recoil energies only the low den- 
sity flows, with high total velocity, can contribute and 
consequently the signal is far smaller than at low recoil 
energies. The lower density flows have negative v$ and, 
in contrast to the high density flows, have larger speeds 
in June than in December, so that at high recoil energies 
the contribution due to a given flow is slightly larger in 
December. As m x is increased the variations in T(E) 
are moved to higher E. The presence of a Maxwellian 
background smoothes the stepped variations in the sig- 
nal produced by the flows, but they are still discernible 
and if detected would provide a distinctive indication of 
the presence of velocity flows. 

In Fig. 3 we plot the differential event rate, dR/dE, 
for the velocity flows plus isothermal background model 
and also for a pure Maxwellian halo, for a Nal detector 
using the best fit values of the WIMP mass and cross- 
section found by the DAMA collaboration, m x — 54 GeV 
£crp = 4 x 1CT 6 pb. 

Brhlik and Roszkowski jl(| have devised a technique 
for comparing the experimental data released by DAMA 
with theoretical predictions for the annual modulation 
signal, in the absence of detailed information about the 
experimental set-up, such as the efficiency of each Nal 
crystal. Their technique, which is effectively a least- 
squares comparison of the experimental data with the 
theoretical predictions, has been used to examine the 
region of mass-cross-section parameter space compat- 
ible with the DAMA results for various simple non- 
standard, but close to Maxwellian, halo models JlO| , p^ |. 
The best fit values and errors for Sq^ and S m> k released 
by DAMA are calculated under the assumption that the 
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FIG. 3. The differential event rate dR/dE in June (solid 
line) and December (dotted line) for the caustic ring halo 
model with velocity flows plus a Maxwellian background 
as described in the text, and for a pure Maxwellian halo 
(June- long dashed line, December-short dashed line), for 
WIMP mass m x = 54 GeV and cross-section £cr p = 4 x 10~ 6 
pb, as found by the DAMA collaboration, for a Nal detector. 

recoil energy spectrum has the shape expected from a 
Maxwellian velocity distribution. Therefore whilst Brhlik 
and Roszkowski's technique can be used for halo models 
which produce recoil energy spectra close to that pro- 
duced by a Maxwellian, it can not be applied to the caus- 
tic flow model. Furthermore Gelmini and Gondolo have 
recently found that at low recoil energies the annual mod- 
ulation produced by the flows is poorly approximated by 
a sinusoidal . 

We can therefore only make a qualitative discussion of 
the effect of the presence of velocity flows on the DAMA 
allowed region. Brhlik and Roszkowski jl(| found that, 
for a pure Maxwellian halo, the cut-off at large WIMP 
masses in the allowed region is determined by the time- 
dependent part of the signal (i.e. S^ k ), whilst the lower 
limit on the WIMP mass depends on both the time inde- 
pendent and dependent parts of the signal. When a ve- 
locity flow component is added to the Maxwellian back- 
ground the recoil energy spectrum falls off less rapidly 
with increasing recoil energy for large WIMP masses, 
whilst for smaller WIMP masses the recoil energy spec- 
trum falls off more rapidly with increasing energy (see 
Figs. 2 and 3). This suggests that the range of WIMP 
masses compatible with the energy distribution observed 
by DAMA would be likely to be smaller for the caustic 
ring model than for a pure Maxwellian halo. The allowed 
region obviously also depends on the magnitude and sign 
of the annual modulation. For the caustic flow model 
the sign of the modulation is opposite to that observed by 
DAMA for some, WIMP mass dependent, ranges of recoil 
energy, however since the experimental data is binned in 
IkeV bins this may not prevent the velocity flow model 
being consistent with the DAMA data. It is possible 
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though that the distinctive effects of the velocity flows on 
the recoil energy spectrum and on the sign of the annual 
modulation could lead to limits on the allowed fraction 
of the local halo density in velocity flows. 

For the purpose of estimating WIMP direct detection 
rates the assumption of a standard Maxwellian halo is 
certainly reasonable. Now that experiments are reach- 
ing the region of parameter space populated by super- 
symmetric models, and in the case of DAMA claiming a 
positive signal, it is important to extend the theoretical 
analysis to more sophisticated, and hopefully more real- 
istic, halo models. This process will be facilitated by the 
public release of data in a form subject to the minimum 
number of theoretical assumptions possible. 



V. CONCLUSIONS 

In this paper we have studied the WIMP direct detec- 
tion signal, in particular its annual modulation, for the 
caustic ring halo model. In this model the WIMP distri- 
bution at the Earth's location has a series of peaks in ve- 
locity space, corresponding to particles which are falling 
into the galaxy for the first time and those which have 
fallen in and out a number of times but have not yet been 
thermalized. These peaks produce a distinctive imprint 
in the differential event rate, with the sign of the annual 
modulation (i.e. whether the event rate is larger in June 
or December) changing with detector recoil energy. The 
presence of an isothermal background component to the 
halo smoothes out the sharp changes in the differential 
event rate produced by the velocity flows but the distinc- 
tive changes in the sign of the annual modulation remain 
potentially discernible. 

Finally we discussed the compatibility of this model 
with the results of the DAMA experiment. The recoil en- 
ergy spectrum varies more rapidly with WIMP mass than 
that produced by a standard Maxwellian halo, whilst for 
some recoil energies the annual modulation signal has 
the opposite sign to that observed by DAMA. These ef- 
fects suggest that for this model the region of WIMP 
mass-cross-section parameter space compatible with the 
DAMA data would be smaller than for the standard 
Maxwellian halo model. In addition it may be possible, 
via a full likelihood analysis, to constrain the fraction of 
the local halo density in velocity flows. This illustrates 
that if a significant component of the galactic dark mat- 
ter is composed of WIMPs, then WIMP direct detection 
experiments with fine grained directional and energy res- 
olution may be able to probe the local galactic structure, 
complementing the information which indirect detection 
experiments |28l would be able provide on larger scales. 
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